###Anlise exploratria de dados

getwd()
setwd("C:/Users/Edmundo/Documents/Meus Documentos/BIANCA/MESTRADO/Linguagem R")
getwd()

##Cervejas

cervejas <- c("chope","lata","garrafa","chope","garrafa","garrafa","lata","lata","nenhuma","lata","garrafa","garrafa","garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma","garrafa","garrafa","garrafa","chope","garrafa","garrafa","garrafa","chope","garrafa","lata","lata")
help(barplot)
help(dotplot)
barplot(table(cervejas))
dotchart(as.matrix(table(cervejas)))
#O que tem maior razo dado/tinta  o grfico de pontos.

##Rios

rivers
data(rivers)
proporcao <- length(rivers[rivers<mean(rivers)])/length(rivers)
proporcao #66.67%
quantile(rivers,prob=0.75) #680
mean(rivers)
mean(rivers,trim=0.25)
median(rivers)
#Existem sim diferenas ente esses valores, a mdia truncada est mais prxima da mediana do que da mdia, representando a existncia de valores extremos na determinao da mdia.

##Caixetais

caixeta <- read.csv("caixeta.csv",as.is=TRUE)
caixeta
str(caixeta)
par(mfrow=c(2,2))
hist(caixeta$h)
hist(caixeta$h[caixeta$local=="chauas"])
hist(caixeta$h[caixeta$local=="jureia"])
hist(caixeta$h[caixeta$local=="retiro"])
#Existe diferena entre a distribuio dos tamanhos dos caixetais de acordo com a localizao. Jureia e Retiro parecem ter uma distribuio menos dominante de rvores de acordo com a altura dos fustes.

##Eucaliptos

eucalipto <- read.table("egrandis.csv",header=TRUE,sep=";",dec=",",as.is=TRUE)
eucalipto
str(eucalipto)
boxplot(as.numeric(eucalipto$dap)~as.factor(eucalipto$regiao))
boxplot(as.numeric(eucalipto$dap)~as.factor(eucalipto$rot))
qqnorm(as.numeric(eucalipto$h))
qqline(as.numeric(eucalipto$h))
#A distribuio dos dados no segue uma curva normal.

##Mais caixetais

par(mfrow=c(1,1))
plot(caixeta$cap[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="retiro"]~caixeta$h[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="retiro"])
plot(caixeta$cap[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="jureia"]~caixeta$h[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="jureia"])
plot(caixeta$cap[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="chauas"]~caixeta$h[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="chauas"])
scatter.smooth(caixeta$cap[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="retiro"]~caixeta$h[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="retiro"])
scatter.smooth(caixeta$cap[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="jureia"]~caixeta$h[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="jureia"])
scatter.smooth(caixeta$cap[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="chauas"]~caixeta$h[caixeta$especie=="Tabebuia cassinoides"&caixeta$local=="chauas"])
library(lattice)
xyplot(caixeta$cap[caixeta$especie=="Tabebuia cassinoides"]~caixeta$h[caixeta$especie=="Tabebuia cassinoides"]|caixeta$local)